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We construct a mean-field formulation of the thermodynamics of ion solvation in immiscible polar 
binary mixtures. Assuming an equilibrium planar interface separating two semi-infinite regions of 
different constant dielectric medium, we study the electrostatic phenomenon of differential adsorp- 
tion of ions at the interface. Using general thermodynamic considerations, we construct the mean- 
field f2-potential and demonstrate the spontaneous formation of an electric double-layer around the 
interface necessarily follow. In our framework, we can also relate both the bulk ion densities in the 
two phases and the distribution potential across the interface to the fundamental Born free energy 
of ion polarization. We further illustrate this selective ion adsorption phenomenon in respective 
examples of fully permeable membranes that are neutral, negative, or positive in charge polarity. 
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I. INTRODUCTION 

Electro-chemical and physical processes at liquid-liquid 
interfaces and across membranes are broadly important 
in many systems. Examples range from cellular and phys- 
iological systems [Tj, Q to everyday applications, such as 
portable batteries, as well as scenarios for the origin of 
life 3 and also biotechnology Q . Biologically, ion trans- 
port is essential for transmembranous and transcellular 
electric potential, fluid transport, and maintaining cellu- 
lar volume [? ] ; the failure in the regulation of the above 
would lead to such conditions as septicemia induced pul- 
monary edema, electrolyte abnormalities in pyelonephri- 
tis of early infancy, hypovolemia and hyponatremia 6]. 
In the chemistry community, there has been a long in- 
terest in the rich electrochemistry associated with inter- 
facial charge transfer that is important in, for example, 
hydometallurgy and two- phase electrolysis Q. 

The interaction of ionic species, particularly pro- 
teins and nucleic acids, with interfaces and membranes 
has become exceedingly iinportant in our understand- 
ing of biological systems [l[ and in processes such as 
genetic transformation, since the latter specifically in- 
volves DNA crossing cellular membranes. Moreover, im- 
portant biomolccular processes take place more efficiently 



•Electronic address: 'w-kung@northweste rn.edu| 

t Electronic address: francisco.solis@asu.edul 

■f Electronic address: im-olvera@northwestern.edul 



at heterogeneous media than in homogeneous systems, 
as demonstrated in the example of DNA renaturation at 
water-phenol interfaces |^]. Therefore, the study of ionic 
profiles at interfaces also leads to a better understand- 
ing of these various biomolecular processes. In general, 
in situ experimental techniques are required to evaluate 
interface phenomena. For example, x-ray-standing waves 
measurements provide theprofile of ions j8|] and also more 
recently of nucleic acids , both along hard surfaces in 
aqueous solutions p^ . On the other hand, recent exper- 
imental techniques have detected ion density profiles at 
liquid-liquid interfaces with the use of x-ray reflectivity 
measurements [ill . [l2j . These experiments provide di- 
rect knowledge of electrostatic interactions among ionic 
species at interfaces and can be easily extended to ana- 
lyze the adsorption of proteins and nucleic acids at liquid- 
liquid interfaces. 

In their experimental study [Til . [T^ , Luo et al. demon- 
strated that the predictions from Gouy-Chapman theory 
did not match well with x-ray reflectivity measurements, 
presumably due to its neglect of molecular-scale struc- 
ture in the liquid solution. To better fit their experi- 
mental data and the free-energy profile of the ions, the 
authors instead formally introduce an ad hoc term in the 
electrostatic energy, which presumably includes effects 
from short-range correlations and which can be theoreti- 
cally approximated by the potential of mean force, expe- 
rienced by single ions near the interface, obtainable from 
molecular dynamics (MD) simulations. Notwithstanding 
the absence of a fundamental theoretical framework for 
the aforementioned system of a planar liquid-liquid in- 



2 



terface, progress has been made with phenomenological 
approaches, including the example of a Ginzburg-Landau 
theory of ions solvation and their distribution around an 
interface [13j . 

In this work, we will present a top-down thermody- 
namic formulation for the system of dissolved ions in im- 
miscible polar-binary mixtures. In our formulation, we 
can naturally incorporate details such as the interaction 
between ions and background solvents, at the mean-field 
level, via the Born approximation [3, UB] • We will then 
relate the various hitherto phenomenological parameters, 
namely, the bulk ion densities in each phase as well as the 
distribution potential across the interface, to the funda- 
mental parameters of our system that include the dielec- 
tric constant of each phase and the size and charge of 
dissolved ions. In addition, the full nonlinear Poisson- 
Boltzmann equation follows from variational principle 
and can be analytically solved for the electrostatic po- 
tential and for the ion density profiles without further 
approximations. Electro-neutrality, in our formulation, 
simply manifests from consistency requirement with ther- 
modynamics. 

As an important consequence of thermodynamics and 
electrostatics, we will demonstrate the physics of differ- 
ential adsorption of ions near an interface, which is a gen- 
eral phenomenon occurring between two phases of differ- 
ent dielectric medium wherein ions are selectively driven 
into confinement near, or expulsion from, the interface 
based on their charge polarity. In particular, we will 
consider examples involving a neutral membrane, a nega- 
tively charged membrane, as well as a positively charged 
one partitioning two chemically different environments. 
Given the long-standing interest in this topic across di- 
verse disciplines and the many still unresolved issues, it is 
inevitable that there exists many prior works on the sub- 
ject with very different emphases and approaches. For 
example, chemists have traditionally focussed on detail- 
ing the microscopic molecular interactions between ions, 
solvents, and the different types of membranes 0,[lB,[i3- 
For the same reason, there exists considerable varied us- 
age of terminology 'la\ regarding similar concepts. For 
clarity, we will henceforth explicitly consider in this work 
only permeable membranes. In addition, to establish no- 
tations and conventions, we will refer to the potential 
gradient across the interface, due to the different exist- 
ing bulk ions densities in the two phases, as the distri- 
bution potential ^d- The cases of neutral and charged 
membranes (regardless of polarity) would correspond to 
the terminology of nonpolarizable and polarizable liquid- 
liquid interfaces, respectively, as discussed in chemistry 
literature. While partial results regarding these systems 
have been derived in various contexts and guises, our 
present work represents a unifying framework in a com- 
plete, mean-field thermodynamic formulation for the sys- 
tem of dissolved ions in immiscible polar binary mixtures, 
within which we can derive rigorous relations between 
the different phenomenological parameters of the system 
that would otherwise be unrelated when treated in the 




FIG. 1: A single charged plate in salt solution. The system is 
symmetric with respect to z — 0. Without loss of generality, 
the charged plate is assumed to have a positive surface-charge 
density a, with counterions (solid circles) and coions (open 
circles) present on both sides of the plate. 

context of electrostatics alone. 

In what follows, we will present, in section II, our 
thermodynamic formalism in terms of a mean-field fl- 
potential and derive, via variational principle, the rel- 
evant electrostatics and thermodynamic constraints in 
our system. We will then proceed in section III to con- 
sider full solutions to the nonlinear Poisson-Boltzmann 
equation in this new context and demonstrate the gen- 
eral mechanism of differential adsorption of ions across a 
liquid-liquid interface. Section IV concludes our paper. 

II. THERMODYNAMICS OF ION SOLVATION 

For the system of dissolved salt ions in polar binary 
mixtures, there are several parameters that have been, 
as previously mentioned, traditionally considered as inde- 
pendent inputs to the formulation of such systems. These 
input parameters include, for each phase a, the dielectric 
constants £„, the bulk ion densities Poj-, and the dis- 
tribution potential across the interface To further 
characterize these parameters, we now turn to thermody- 
namic considerations. We will henceforth denote the two 
polar phases by a = A, B with A labeling the more polar 
phase. In principle, the condition of equilibrium implies 
equal chemical potential across the interface, /i^ = /i^, 
which should in turn determine the bulk densities of the 
salt ions existing in each of the two phases. The thermo- 
dynamics with this equilibrium condition are formulated 
by constructing the fi-potential (equals to —PV where 
P is the pressure and V is the volume) for each phase a 
as follows: 

= + + (i) 

where we have divided the potential into a contribution 
due to electrostatics fiej, to thermodynamics ^th, and to 
solvation interactions Qgoi- In what follows, all energy 
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quantities will be measured in units of fc^T, all length 
scales in units of the Bjerrum length — /ksT, and 
all ionic charges quoted as multiples of the fundamental 
electron charge e. The electrostatic part delineates the 
free-energy dependence on the dielectric constant e" and 
imposes Gauss' law via a Lagrange-multiplier constraint 
using the electrostatic potential $: 



' (2) 

On the other hand, the thermodynamic contribution con- 
sists of the entropy associated with the ions, treated as 
point particles, in solution as well as their corresponding 
chemical potentials /i^ , 



the various factors within the integral, it is straightfor- 
ward to see that r2";|fc=o V^^'^, for poc 7^ 0, which is 
inconsistent with the thermodynamic construction of the 
extensive fi-potential. Therefore, the requirement of con- 
sistency with equilibrium thermodynamics automatically 
imposes charge neutrality in the overall system; namely, 
Pq = (we note that this elelctroneutrality condition in 
linearized approaches appears as a restriction to avoid 
the divergence in the potential (lOj). 

From Eq. ([T]), variation with respect to the electric 
field E"-, the Lagrange multiplier (electrostatic poten- 
tial) and the number density readily yields the 
relations 



(5^;" : £;°(a;) = -V$"(£c), 
"e"£;"(a;) 



(6) 



(5$ 



V • 



An 



p^+(a^)+p^(a^), (7) 



\ogpl{x)±^^{x)-gl^^,l. 



(8) 



17^/. = P+(a;)log 



p%(a 



P- {x) log 



pt{x) 



-p%p%{x) 



piplix), 



(3) 



where the chemical potentials will be specified by the 
equilibrium constraint that they should be equal across 
the interface. Lastly, the solvation interaction between 
the ions and the respective polar solvent in each phase 
is modeled by the mean-field Born approximation of po- 
larziation; 



where the Born polarization energy is given by 



(4) 



We note that Eqs. ([6|) and ([7| simply reproduce, respec- 
tively, the expressibility of the electric field E in terms of 
a scalar potential <i>" and the Gauss' law in electrostatics. 
The last relation in Eq. ([8]) provides an explicit expres- 
sion of the chemical potential for each phase in terms 
of its respective input parameters such as the dielectric 
constants and ion sizes. It is now obvious that the ther- 
modynamic condition of equal chemical potentials across 
the different phases a would provide further constraints 
on the respective bulk ionic densities. To make connec- 
tion with the distribution potential we also need to 
introduce explicitly now the further assumption of Boltz- 
mann distribution for the ion densities. In particular, we 
have, for the more polar phase A, 



9± 



(5) 



for ions with valency Z± and radius R±. As expected, 
we note that an increase in dielectric constant e" or in 
the ionic radius R± for a given valence Z± facilitates sol- 
vation of ions in polar mixtures and leads to an overall 
lower Born energy. Consequently, Born's model predicts 
a larger solubility for anions (which are generally larger in 
size due to reduction in redox chemistry) than for the cor- 
responding cations. As we will see in subsequent sections, 
both charge properties and the size factor play an impor- 
tant role in the physics of ion solvation near a boundary 
interface. 

To examine the consequence of our thermodynamic 
formulation of ion solvation in immiscible polar binary 
mixtures, we first perform dimensional analysis on the 
electrostatic term. We note that the fi-potential is an 
extensive quantity by construction: oc V, where 

V is the volume of the overall system. Defining the 
total charge density p"{x) = |Z+|p"(a;) — \Z^\p'^{x), 
we assume that it has a nonvanishing zero mode and 
observe that its contribution to the energy density is 
{Pocf S J dy i ^^^i . Based on the dimensionality of 



pHx) 



Po e 



(9) 



while a similar expression (modulus the exponential fac- 
tor of the distribution potential $d) applies to the less 
polar phase B. Taking the logarithm of Eq. ^ readily 
yields 



logp^(cc) ±$^(a;) = logp^±$D 



(10) 



Again, a similar expression to Eq. (jlOp . modulus the 
term containing the distribution potential, applies to the 
less polar phase B. Now upon substituting Eq. (fTn|) 
and the corresponding expression for phase B into Eq. 
([5]), as well as imposing equal chemical potentials across 
boundary, we obtain the following relations 



^ogPo 



A ^_ „A _ i„„ B „B 



B 



\ogp{^ ~^D- g- = \og p" -gi 



(11) 
(12) 



We can now take the sum and difference of the two equa- 
tions to obtain the relations 
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.93, 



(13) 
(14) 
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FIG. 2: Comparison of the electrostatic potential across 
a positively and negatively charged membrane in nitroben- 
zene/water solution at 68F. The surface charge densities are 
quoted in units of tAi^A and $^ ~ 0.112. The potential 
is plotted for a = -1.5(A), a = -0.93(n), a = -0.5(V), 
C7 = 0(o), C7 = 0.15(A), cr = 0.32(B), a = 0.5(T). 



where we have defined the average Born polarization en- 
ergy over ions of both charge polarity, ga = ga+ + <7q-, 
and the corresponding differential Born energy, A^q, = 
(7a+ — <?a- between the cations and anions. In terms of 
these two quantities, we can relate the bulk ion densities 
and the distribution potential to the dielectric constants 
£a in each phase via 



P± 



(15) 
(16) 



III. ELECTROSTATICS OF DIFFERENTIAL 
ADSORPTION 

As mentioned in the Introduction, it is now feasible 
to make precise experimental detection of the ion dis- 
tributions near liquid-liquid interfaces [111, [l2|. It has 
been found that while the mean-field Poisson-Boltzmann 
framework provides a generally good description for the 
ion profiles that matches well with experimental re- 
sults obtained by x-ray structural measurements, the lin- 
earized Guoy-Chapman approximation is insufficient for 
most general cases [l2|. For completeness, we will now il- 
lustrate the solution to the Poisson-Boltzmann equation 
for the case when the membrane across the two dielectric 
media is charged (polarizable), including the particular 
interesting case of the neutral (non-polarizable) interface. 
Furthermore, our liquid-liquid interface is treated as an 
infinite two-dimensional plane so that we do not need to 
consider surface terms. We thus essentially reduce our 
system to a one-dimensional problem. Denoting the lat- 
eral distance from the membrane by z, we will now for- 
mulate the full non-linear Poisson-Boltzmann equation 
with the proper boundary conditions and obtain the cor- 
responding electrostatic potential $"(z) and ion density 
profiles p±{z) in each phase a. 

To begin, we analytically formulate our system in 
which the two phases A and B are in equilibrium, where 
each phase is characterized by its respective constant di- 
electric permittivity ea and sb, as follows: 



eaV-E^ = P^, z>0, 
sbV-E'^ - pf, z<0, 



(17) 
(18) 



subject to the following boundary conditions at z = 0: 



where A^Afji — Ag^ — Ag^- 

In summary, we have thus shown in Eq. (|15p that the 
bulk densities Pq_|_ are related to the difference of the 
average Born solvation energy between the two phases, 
while in Eq. ([TBI) the distribution potential is related 
to the difference in the differential Born energy between 
ion types in the two phases. We note that even though 
the above formulation has been done for a neutral liquid- 
liquid interface wherein the charge distribution of the 
system resides completely amongst the dissolved ions in 
the bulk, it is straightforward to generalize to the case 
where a surface charge exists on the interface itself. 

To determine the equilibrium ion density profiles {x) 
as well as the electrostatic potential $(a;), it remains for 
us to solve the Poisson equation [Eq. ([7])] after apply- 
ing the intcgrability condition [Eq. ©] as well as the 
thermodynamic ansatz of the Boltzmann distribution for 
the ion density [Eq. ([9])]. We will illustrate with two 
classes of examples in the subsequent section involving 
non-polarizable (neutral) but permeable membranes as 
well as polarizable (charged) membranes. 





= Ana, 
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= Eyl^^Q- . 


(21) 



Making use of Eq. ^ and substituting the Boltzmann 
distribution for ion density profile [Eq. ([9|)] , we now take 
advantage of the one-dimensional nature of our system 
and collapse the above formulation into a single equation, 
written in terms of the electrostatic potential $"(z), to 
obtain the Poisson Boltzmann-equation representation of 
the Gauss' law: 

^ - (K^)%inh($^-<i>D) , z<0, (22) 
= (K^)'sinh$^, z>0, (23) 

where, respectively, — ^ Aire^p'^ / e'^kT is the inverse 
Debye screening length, and Pq = Z\pq + Z'^p^ _ is 
the combined bulk densities of the ions in the two phases. 




FIG. 3: Comparison of the ion distributions of ions in the two polar phases, A{z < 0) and B {z > 0) characterized by ea and 
eb, respectively, when the surface density of the interface is cr = —1.5(A), a = — 0.93(n), a — — 0.5(V), a = 0(o), a — 0.15(a), 
a = 0.32(B), a — 0.5(T). (1) Cations. (2) Anions. (3) Total charge density p = p+ + p-. The curves are presented in separate 
plots because of difference in concentrations. For this particular plot, we use a range of surface charge densities a in units of 
EA and — 0.112 for the nitrobenzene/water solution at 68F. 



a = A, B. The inverse Debye screening length has been 
so defined as to scale the electrostatic potential di- 
mensionless. The corresponding boundary conditions for 
the potential are as follows: 
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(24) 


dz 


z^o- dz 


— Aira, 


(25) 




d^^ 

dz 


z — >~oo 


= 0, 


(26) 




d'P^ 

dz 


z — ^+oo 


= 0. 


(27) 



The general solution to the full nonlinear Poisson- 
Boltzmann equation can be readily found [l^l. At 
present, we simply furnish the particular solution for our 
system of nonpolarizablc liquid-liquid interface specified 
by the above boundary conditions: 

$^(z) = -4tanh"^ (C^e'^-^^) , z < 0, (28) 



$^(z) = 4tanh"^ {Cse- 



As shown, Eqs. 



and 



2 > 0. (29) 



are expressed in terms 



of the two integration constants, Ca and Cb, and as 
such, physically meaningful solutions would correspond 
to the range of \Ca\ < 1. Alternatively, we can rewrite 
Ca — exp(— KqZq), and the solutions now take the 
form of ^^{z) = - 4sgn(CA) tanh""^ ^ka(z-za)^ g^^^ 
^^(z) = 4sgn(C_B) tanh^"^ e~'*^(^~"'^\ where sgn[ • ] is 
the sign function. The choice of signs in the above solu- 
tions is such that the integration constants are positive 
when we have a neutral interface, a = 0. 

In terms of the integration constants Ca and Cb, the 
continuity equation in Eq. (|24p reads: 



(I + CaKI + Cb) 



e 2 



(1 - Ca){1 - Cb) 
while the Gauss' law in Eq. (P5)l takes the form of 
AcakACa AeBKBCB 



(-2 



1 



f<2 



1 



— Ana. 



(30) 



(31) 



This system of equations admits the following simple al- 
gebraic solutions for Ca and Cb- 

EAKA + eBKB - R 
= _ , — ^ — 77; • (o2) 



Ce 



27ra + €bHb sinh $15/2 
^B^B + ^aka cosh$£)/2 — R 
— 27r(T -|- e^iKA sinh $13/2 ' 



(33) 
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where 



+ 2eAKAi^BKB C0sh$£)/2. 



(34) 



Thus, we have now determined the constants of integra- 
tion expHcitly in terms of the input parameters ca, es, 
and As such, we can write down the positive and 
negative ion densities in a straightforward manner: 



A l'l±CAe^ 
Po 



Po 



1 T Cab'^^ 
1 ± Cse-^s; 



(35) 
(36) 



To iUustrate the profiles of the electrostatic potential 
$(2;) and of the ion distributions p%(z) across an inter- 
face, we will now consider, in particular, the system of 
a polarizable membrane between a nitrobenzene solution 
= 35.7£o at 68F] of sodium chloride and a water so- 
lution [e^ = 80.4eo at 68F] of the same salt. Given the 
size ratio between the ions of k. 0.695, we ob- 

tain the following normalized relations between the var- 
ious Born solvation energies: g'^ w 0.662, « 0.446, 
1.44, and g^ = 1. Thus, it follows readily that 



9l 



p^IpE 



,1.33 



« 3.79 and <I>d k, 0.112, which implies 
1.30. We display the resulting plots in 



that k^/k^ f 
Figs, m and O 

The excess charge, Q%, defined as 



(37) 



and due physically to the accumulation or depletion of 
charges near the interface as the ions migrate across the 
boundary, can also be obtained in analytical closed-form 
for the present case of planar interface. We note that the 
excess charge (5± has the dimensions of charge density 
per unit of transverse area. Upon direct evaluation, we 
obtain 

1 eA^ACA 

27r 1 T ' 
1 cbkbCb 



Qi 



27r \±Cb 



(38) 
(39) 



The total excess charge in each phase a can be simply 
evaluated by = -I- . In addition, the net charge 
flux across the interface can be defined as AQ = ^{Q^ — 
Q^). The total net charge accumulated in each phase 
compensates for the charge introduced by the polarizable 
membrane to the overall system such that the relation 
CT + Q^ + Q^ ={) holds at all times. The plot of Q^, 
and AQ for various values of the surface charge density 
a is shown in Fig. [H 

In the case of polarizable membranes, we see that the 
additional surface charge density can have dramatic ef- 
fects on the electrostatic potential and charge distribu- 
tions of the two phases. As shown in Fig. [51 there exists a 
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FIG. 4: Plots of the total excess charge , and the net 
charge flux across the interface AQ. For this particular plot, 
we use a range of surface charge densities a in units of £a ka 
and = 0.112 for the nitrobenzene/water solution at 68F. 



sufficiently large surface charge density on the interface, 
cr*", at each phase a such that the electric field vanishes 
as manifested by the corresponding constant electrostatic 
potential. The further addition of surface charges, either 
positive or negative, reverses the sign of the potential 
gradient and of the excess charge in the corresponding 
phase. The critical surface charge densities are found to 
be 



1 . 

—eBKB smh-— , 

ZTT Z 
1 . , $13 

-TT'^A'^A smh—. 

ZTT Z 



(40) 
(41) 



It can be also checked that at these critical surface charge 
densities, our expressions for the excess charge Q" and 
the net charge flux across the interface AQ also vanish in 
the respective phase (Fig. 4). Notwithstanding the seem- 
ing linearity that both Q" and AQ follow with respect to 
the surface density a, we do observe a slight jump across 
(7 = thereby establishing a small deviation to perfect 
linearity. 

We now consider a few general properties of our solu- 
tions. Firstly, we note that the distribution potential 
across the interface, itself generated by the different di- 
electric constants of each phase, is only dependent upon 
ion types and not of their concentration. Thus, in the 
absence of any ions in the system, we must then recover 
a uniform potential. Indeed, we observe exactly that in 
approaching the limit of zero ionic concentration where 
the gradient of the potential across the interface gradu- 
ally decreases; the required net jump is being pushed 
outwardly towards infinity. 

Another noteworthy feature is the asymptotic behavior 
of the potential at large distances from the interface. In 
the limit of z ^ 00, we have <I>^ « — AC a exp {ka z), 
and $^ « 4Cb exp (— kbz). Thus, we see that the in- 
tegration constants now act as potential sources in the 
solution form for the linearized version of the Poisson- 
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Boltzmann equation. Given that each integration con- 
stant satisfies the constraint of |Cq| < 1, the asymptotics 
always arises from a hnear source of maximum potential 

^max ^ 4. 

Lastly, we remark that our results are consistent with 
the Born description of solvation energy: since cations 
are in general smaller in size due to their oxidized state, 
the solvation free energy is always more negative when 
compared with their anionic counterpart, provided that 
they are both subjected to the same solvent of the same 
constant dielectric permittivity. Therefore, it is always 
thermodynamically more favorable to dissolve more of 
the cations in the more polar phase and redistribute the 
density profiles of the remaining ions accordingly in both 
phases. The validity of our argument remains even when 
the membrane becomes charged in the same polarity as 
the cations. In general, we do not expect that the multi- 
valency of the dissolved ions would qualitatively change 
the behavior of the differential adsorption phenomenon 
described herein, at the mean-field level, other than to 
intensify its effect as evident in the prefactor of our Born 
expression of solvation energy. 

IV. CONCLUSIONS 

In this work, we have constructed a self-contained ther- 
modynamic formulation of ion solvation in binary im- 
miscible polar mixtures and described the general phe- 
nomenon of differential ion adsorption at liquid-liquid in- 
terfaces via the full nonlinear Poisson-Boltzmann frame- 
work. Assuming the Born model of solvation and us- 
ing the mechanics of the fi-potential, our formulation 
presents a fundamental, mean-field description that re- 
lates the experimentally detectable bulk-ion densities in 
each phase, as well as the distribution potential across 
the membrane, to the respective dielectric constants and 
sizes of each ion species present in the two phases. We 
note that electrostatics alone is insufficient in capturing 
the physics of ion solvation and interfacial adsorption. 

Our work can be generalized in several ways. Our for- 
mulation models the dissolved ions specifically as point 
particles and the liquid-liquid interface as sharp. Im- 
miscible multicomponcnt liquids in some cases can lead 
to the broadening of interfaces [2i|. In such instances, 
adsorption of ionic component in slabs of immiscible liq- 
uids has been recently considered using liquid-liquid weak 
segregation approximations pLSj . Thus, it is of interest to 
consider other morphologies and systems than those pre- 
sented in this work. As a refinement to our model, it is 
also possible to incorporate further molecular details in 



the construction of our thermodynamic potential. This 
can be done, for example, with the explicit inclusion of 
short-range interactions between the ions by replacing 
our expression of f2"^; by 

^Toi = J J drdr'Gsr{r-r')p{r')p{r) (42) 

One could then consider components with more com- 
plicated structures than point ions, such as polyelec- 
trolytes, and consider the additional effects introduced 
by these new degrees of freedom [l^, [2^ . Though pro- 
tein adsorption on hard surfaces have been analyzed [24| , 
a general description of macroion adsorption to liquid 
interfaces is lacking. Such a general description would 
be highly useful, as, for example, the phenomenon of 
DNA adsorption to the interface among liquid surfaces 
is crucial in cell biology [l|, Q • The inclusion of proteins 
and/or DNA molecules in the system and the study of 
their corresponding adsorption phenomenon along liquid- 
liquid interfaces will be presented elsewhere. We point 
out in passing that the long-range interaction of electro- 
statics would not contribute at the mean-field level due 
to overall electro-neutrality. 

It is our hope to have demonstrated in this work 
that the differential adsorption of ions is a general phe- 
nomenon in electrostatics that occurs along interfaces be- 
tween different dielectric media. Such chemical environ- 
ments are ubiquitous in many biological systems at the 
cellular and physiological levels In particular, it was 
pointed out by Onsager Q regarding the potential impor- 
tance of such a "primodial oil slick" in the early devel- 
opment of life. The phenomenon of thermodynamically 
driven, selective confinement of ions and the consequent 
breaking of charge-conjugation symmetry near a liquid- 
liquid interface would provide just the physical mecha- 
nisms necessary for the essential organic molecules to ag- 
gregate in close, two-dimensional proximity to each other, 
where the process of diffusion would have worked much 
more efficiently in bringing together these molecules and 
starting the chain reaction of life than in the otherwise 
three-dimensional scenario in the bulk. It is our further 
hope to generalize our results in this work to other in- 
terfacial geometries such as the cylindrical and spherical 
cases. 
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